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Abstract 

After summarizing three versions of trajectory-based quantum me¬ 
chanics, it is argued that only the original formulation due to Bohm, 
which uses the Schrodinger wave function to guide the particles, can 
be readily extended to particles with spin. To extend the two wave 
function-free formulations, it is argued that necessarily particle tra¬ 
jectories not only determine location, but also spin. Since spin values 
are discrete, it is natural to revert to a variation of Bohm’s pilot wave 
formulation due originally to Bell. It is shown that within this for¬ 
mulation with stochastic quantum trajectories, a wave function free 
formulation can be obtained. 


1 Introduction 

In the last ten years or so, there has been a growing interest in Bohm’s 
trajectory based interpretation of quantum mechanics PE], Even 
though the introduction of classical-like particle trajectories into the 
quantum mechanical state description does not have any observable 
consequence, it offers more practical ways to depict and explore quan¬ 
tum behavior. Specifically, in quantum chemistry this has led to novel 
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ways to numerically solve and investigate multi-particle quantum sys¬ 
tems [2]. From a conceptual or ontological point of view, Bohrn’s 
interpretation has a number of well-known advantages; the most com¬ 
pelling feature being the resolution of the measurement problem: un¬ 
like the traditional Copenhagen interpretation, there is no need for 
a collapse of the wave function. Obviously Bohm’s trajectory based 
interpretation is not generally adopted and one can raise various ob¬ 
jections and queries. For example: Why is there still a wave function 
in this picture - why should the particles need to be guided by a pi¬ 
lot wave? Why would there only be trajectories for position (or for 
fermion number, if we adopt Bell’s version of Bohm’s interpretation 
i)? The experimental verification of Bell’s inequalities j4] has ef¬ 
fectively excluded local hidden variable interpretations of quantum 
mechanics. The presence of the pilot wave and its manifestation as a 
quantum potential are presumed to provide the required non-locality 
that allows compliance to Bell’s inequalities, but how does this come 
about? 

The first question has been addressed by Holland, Poirier and Hall 
et al. in recent work mum- Interestingly, these authors show that 
Bohm’s quantum trajectories can be obtained without a guiding wave 
function from a reformulated theory that prescribes the equations of 
motion for the particles, along with a probability distribution for the 
resulting particle trajectories. This probability distribution (or, in 
the version of ref. [7j, the repulsion among trajectory realizations) 
generates an additional force that causes all required quantum effects. 

The second question was addressed in [8], where it was shown that 
it is in fact possible to compute trajectories for all observables, includ¬ 
ing inherently discrete entities like spin. In this work it was argued 
that one could choose a preferred complete set of (commuting) observ¬ 
ables, or - lacking compelling arguments to prefer one set over another 
- allow “all” observables, also if they are mutually non-commuting, to 
have well-defined trajectories. 

The third question of how non-locality manifests itself in Bohm’s 
quantum trajectory interpretation, has been discussed in detail in 
ref. pTj . Here it is shown how coherence in the multi-particle wave 
function acts as a non-local guidance for the particle trajectories and 
thus manages to avoid the constraints imposed by Bell’s inequalities. 
Note that the pilot wave, which carries the entangled spin state, is 
essential to make the particles behave such that results of normal 
quantum mechanics are reproduced. 
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This paper revisits the above three points to further clarify the 
issues and show their interconnectedness. First, in section [2] three 
available formulations of Bohm’s trajectory interpretation of quantum 
mechanics for spin-zero particles are summarized. Since all three are 
designed to reproduce observable results of normal quantum mechan¬ 
ics, their difference is primarily in the underlying ontology, or in what 
are assumed to be “elements of reality”. Additional differences, for ex¬ 
ample how amenable each formulation is to numerical evaluation, will 
not be pursued in any detail here. The first formulation, supported 
by Bohm [Tj, Bell [3j and others |9] firmly assumes a pilot-wave, which 
is a solution of the Schrodinger equation, as part of reality; the two 
other formulations make no explicit reference to a wave function or 
Schrodinger equation mmm- 

Then, in section 02 it will be argued that such wave function-free 
formulations of quantum mechanics that only give a “beable” status 
to position, when carried over to particles with spin, will lack the 
state information that is required to provide the quantum correla¬ 
tions required to describe particles with spin. This requires that also 
spin state has to be included somehow as an element of reality, as 
is the case for the stochastic trajectory formulation of quantum me¬ 
chanics developed in ref. [8]. Hence, it suggests to explore if the wave 
function-free formulations discussed in section [2] can be extended to 
this stochastic trajectory formulation. This is carried out in section [U 
where it is shown that the ensemble of quantum trajectories for any 
(discrete) observable can be generated self consistently from a suitably 
chosen distribution of initial values by applying a local (stochastic) 
evolution rule without reference to a wave function. As in the wave 
function-free particle formulation, the dynamics of the stochastically 
evolving quantum numbers is affected by the local potential and non¬ 
local probability distribution defined by the ensemble of trajectories in 
the high-dimensional state space. Finally, some concluding comments 
are presented in Section [6l 


2 Bohm Trajectories without Wave func¬ 
tion: Three formulations 

To keep notations simple and following refs. 13 El, this section will 
consider quantum mechanics of a single particle without spin in one 
dimension. It is mostly straightforward to extend the discussion to 
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non-relativistic quantum mechanics for multiple spin zero particles in 
three space dimensions. The Schroodinger equation for such a one- 
particle system is given by 

ih-ipix, t ) = -^—ip"(x, t) + V(x)ip(x, t), (1) 

2m 

where rn is the particle mass, V the classical potential and Ti Planck’s 
constant. As explained in e.g. [T], for a quantum system described 
by the wave function ijj(x,t) one can define an ensemble of particle 
trajectories x(t,x o). A family of trajectories can be derived from the 
wave function 

$ = Re iS/h . (2) 

using the “pilot” equation, 

x(t,x 0 ) = S'(x 0 ,t)/m, (3) 

were S'(xo,t) is the gradient of the phase of the wave function, eval¬ 
uated at x(to) = xo- It follows from the Schrodinger equation that 
these particle trajectories obey almost classical equations of motion, 

x = ~(V(x) + Q(x,t))'/m. (4) 


Quantum effects are due to the quantum potential that is present in 
addition to the classical potential V, 


Q 


-ti 2 R" 
2m R 


(5) 


Quantum uncertainty arises because one must consider the full ensem¬ 
ble of trajectories, with all possible initial conditions xq, which must 
be distributed initially according to the probability density at start¬ 
ing time to, P(xo,to)- The probability density for particle positions is 
computed from the radial part of the wave function as, 


P(x, t ) = R 2 (x, t). 


( 6 ) 


When evaluating probabilities as averages over this evolving ensemble 
of particles, the same results as in normal quantum mechanics are 
recovered. 

To solve the dynamics of this ensemble of particles, one can follow 
at least three approaches, which will be summarized next. 
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2.1 Trajectories Guided by a Wave Function 
(F-I) 

The most straightforward (and oldest [I] [9] ) approach is to first solve 
the (linear) Schrodinger equation to obtain the complex valued t). 
After computing ip(x,t), one can use its phase S(x,t ) with eq. Q to 
compute particle trajectories. The radial part of the wave function 
does not play a direct role in this formulation, other than providing 
the probability distribution for the initial trajectory positions .To0 


2.2 Trajectories with Self-contained Dynam¬ 
ics (F-II) 


Second, one can exploit the particle aspect of the dynamics more ex¬ 
plicitly and compute particle trajectories from the equation of motion 
([!]) for the particles. Unfortunately, this equation contains the quan¬ 
tum potential, which is only given at the initial time and subsequently 
evolves dynamically as well. In contrast with the previous formulation 
F-I, the quantum potential in this formulation is computed as an emer¬ 
gent property of the ensemble of particles; The function R(x, t ) needed 
to compute the quantum potential © is obtained as the square root of 
the probability density of particle positions in the evolving ensemble 
of particle trajectories. 

Hence, one can apply the following leap-frog scheme for solving the 
particle dynamics: Choose a large set of initial positions, according to 
a (given) probability density, P(xo,to), along with initial values for the 
particle velocities, v(xo,to). Compute the quantum potential at this 
initial time, using ([6]) and ([5]). Make sure that the initial velocities 
satisfy suitable integrability conditions tm, such that they can be 
obtained as the gradient of an underlying scalar function, 

v(x,t 0 ) = S'(x,t 0 )/m. (7) 


Note that this is not straightforward, nor particularly natural, in more 
than one dimension. Then evolve the particle positions in the ensemble 
from to to t = to + dt, using the almost classical equations of motion, 
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h 2 (P 1 / 2 ^,*))' 


v(x,t + dt) = - (y{x)~ 

m 2m P l l z {x,t) 


Ydt , 


( 8 ) 


1 It is even argued mi that the proper probability distribution, P(x,t) = R 2 (x,t) 

will develop dynamically from an arbitrary initial distribution, since only this specific 
distribution is conserved during the non-linear dynamics of the particles. 
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x(x, t + dt) = x(x, t) + (v(x, t + dt))dt. (9) 

Here, the quantum potential is computed from the probability density 
for particle positions, using eq. ((6]). Notice that the velocity is updated 
by adding a gradient, hence condition dTJ) remains valid also for later 
times. After propagating all particle positions in the ensemble to the 
incremented time, the new probability density can be computed using 
a suitable histogramming or averaging technique, 

P(x,t + dt) oc (^particles £ [x,x + a])/a. (10) 

It can be shown m 0 m m that this formulation generates fields 
R(x,t) and S(x,t), which can be computed from P(x,t) and v(x,t), 
and which can be combined into a wave function ijj(x, t ) = R(x, t)e lS ^ x,t ^ h 
that solves the Schrodinger equation ([!]). Hence, as in F-I, this for¬ 
mulation produces the same observable results as normal quantum 
mechanics. This is the case when an infinitely large ensemble of tra¬ 
jectories is used to compute the probability density m- When a 
finite ensemble is used to estimate the distribution, as proposed in [7], 
the quantum potential featuring in (j8]) will be an approximation as 
well, and deviations from normal quantum mechanics will be implied. 

Even though this scheme looks attractive from a numerical point 
of view, because it is relatively simple to solve the time evolution of 
the particles, the obvious drawback is the difficulty in computing the 
time evolution of the quantum potential (or, equivalently, the prob¬ 
ability density of particle positions). One would need either clever 
approximation methods, or a very large ensemble of particle trajecto¬ 
ries - assuming that only very small observable deviations from normal 
quantum mechanics can be allowed. In particular when the number 
of particles increases, the high dimensionality of the state space over 
which the probability density must be evaluated will require a very 
large ensemble size that will make any sampling method quickly pro¬ 
hibitively expensive. 

2.3 Trajectories with Self-consistent Ensem¬ 
ble Dynamics (F-III) 

The third formulation pays tribute to the quantum fluid approach 
originally introduced by Madelung m and more recently picked up 
by the quantum chemistry community mm- It uses the underlying 
particle dynamics, as in section 12.21 above, but use a different method 


6 


to compute the evolving quantum potential. The formulation is best 
explained in discretized from, where the underlying space is replaced 
by a grid. In one dimension, the particle positions are then given by 

x k = ka, (11) 

where a denotes the grid spacing, which for simplicity is taken to be 
constant throughout space and k is the grid cell index. Similarly as 
before, the initial state of the system is defined by specifying for all 
k, the probability density P k (to) for the particle to be located in cell 
k, along with the initial value of the velocity v k (to) that the particle 
would have in this cell. Hence, the state of the single particle is fully 
defined by the values of Pk and Vk for all k. 

As in formulation F-II, the initial values of the velocities in each cell 
must be chosen subject to (the discretized version of) the constraint 
0, and the velocities are updated using the almost classical equations 
of motion 0. The probability density is updated by enforcing the 
continuity equation, 

P k {t + dt) = P k (t) + (P k -i(t)vk-i(t) - P k (t)v k (t))dt/a , (12) 

where it is assumed that the particles can only move to their direct 
neighboring cell, and that the velocity in cell k is in the positive di¬ 
rection. 

In this “quantum hydrodynamics” formulation as it is presented 
above, there is no need to first solve the Schrodinger equation, in 
order to compute the ensemble of particle trajectories. Hence, also 
this third scheme describes particle dynamics without a guiding wave 
function. However, it is also clear that the computational effort to 
simultaneously solve P k and v k for all k could be grosso-modo the 
same as solving the the Schrodinger equation for R and S. 

2.4 Similarities and Differences 

Before closing this section, it may be worthwhile to summarize sim¬ 
ilarities and differences between the three formulations above. It is 
assumed here that the ensemble in formulation F-II (section 12.211 is 
sufficiently large, such that the probability density can be accurately 
computed, and that the grid spacing a in formulation F-III (section 
EH is sufficiently small to avoid discretization effects. For a given 
initial state, all three methods then produce the same trajectories for 
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an ensemble of particles and reproduce all observable results of nor¬ 
mal quantum mechanics. Also, as required 0 ), all formulations have 
non-locality built into their dynamics: in all formulations, particle 
trajectories are influenced by the same quantum potential, which is 
an inseparable function over the multi-particle state space and hence 
intrinsically non-local. 

There will, of course, be differences in computational strategies 
that are most natural or most efficient when implementing a numerical 
scheme for each of these formulations. However, this is not further 
pursued here. There are also important differences in the ontology 
implied by each approach, which will be outlined next. 

In the traditional Bohm interpretation F-I, the wave function along 
with a particle trajectory, are elements of reality. The express pur¬ 
pose of the wave function is to guide the particle along its quantum 
trajectory. Quantum effects, such as interference effects but also en¬ 
tanglement of (spin) states and non-locality, are carried over from the 
wave function, which is evaluated over full state space, to the quantum 
potential [I] or transition probabilities |3j, 

In formulation F-II, there is a very large, possibly infinite, en¬ 
semble of particle trajectories, but no wave function. Particles move 
according to their equations of motion, as in classical dynamics, but 
the equation of motion are endowed with an extra potential (or force) 
term. This (again inseparable) quantum potential is generated as a 
collective effect of the entire ensemble of coexisting state space tra¬ 
jectories. Hence, all trajectories are equally important elements of 
reality. 

For formulation F-III one can adopt an ontology that is close to 
either of the two above, depending on how the ensemble of trajectories 
is interpreted. One can represent the ensemble using the probability 
density and velocity field and give these a similar status as (the mod¬ 
ulus and phase of) the wave function in the Bohm interpretation F-I. 
These evolving probability and velocity fields produce a quantum po¬ 
tential that, along with the classical potential, determines the (single) 
trajectory for the particle. Hence, almost as in F-I, the particle trajec¬ 
tory along with the probability density and velocity field are elements 
of reality. This flavor will be referred to as F-IIIa. Alternatively, one 
can represent the ensemble using a (very large but not necessarily infi¬ 
nite) ensemble of trajectories and effectively adopt the same ontology 
or elements of reality as in F-II. This alternative will be referred to as 
F-IIIb. 


3 Particles with Spin 

Until now, the particles were simple point-particles, with only ‘posi¬ 
tion’ as intrinsic property that is an element of reality. One can argue 
PQ, that it is in fact sufficient that a particle only has a well-defined 
position. Other observable properties, such as momentum or spin, 
need not be intrinsic properties of the particle (and hence elements of 
reality) since they are only indirectly observed though various types of 
measurements that in the end boil down to registrations of pointer po¬ 
sitions or other (displayed or printed) position-type indicators. Even 
though a well-defined spin is not an actual or intrinsic property of the 
moving particle, its trajectory can still be influenced by spin, because 
the guiding wave function evolves differently depending on its spin 
attributes, which are driven by spin-dependent terms in the Hamil¬ 
tonian. In this way the results of the measurement process involving 
particles spin, including particles with entangled states, can reproduce 
the normal results of quantum mechanics [l] [ch. 10,12],[3j [Ch. 18]. 

Since wave function-free formulations require that the dynamics of 
the evolving state can be computed self-consistently without recourse 
to a guiding wave function, this seems to preclude extending formu¬ 
lations F-II and F-III to quantum mechanics for particles with spin, 
unless the particles are also endowed with spin degrees of freedom. As 
Bohm noted |T][Ch.l0.4], this appears to be difficult, since spin de¬ 
grees of freedom linked to a particle would grow proportionally with 
the number of particles, whereas the actual number of spin degrees of 
freedom grows exponentially. Hence, Bohm conceded that for particles 
with spin, a formulation in which trajectories are determined by the 
influence of a combination of classical and quantum forces could no 
longer be maintained If it is possible to avoid Bohm’s conclusion and 
find a wave function-free formulation for trajectories for particles with 
spin, it appears inevitable to start from an ontology that includes full 
multi-particle spin state as an element or reality. Such a formulation 
was presented in ref. [8], which showed how “de Broglie-Bohm-Bell” 
(BBB) trajectories could be computed for any observable - also for 
multi-particle systems with both position and spin. The next section 
will show how this approach can be reformulated to provide a wave 
function-free description of the dynamics. 
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4 BBB Trajectories without Wave Func¬ 
tion 


Following ref. [8], all observables are assumed to take discrete values in 
this section. The Schrodinger wave equation for any discrete quantum 
system, including multi-particle systems with with spin can be written 
as 

iti'ipn — 'y ( H nm 'il > rn , ( 13 ) 

m 

where the composite index n labels the discrete values of the relevant 
observables. For example, eq. (fT3l) can represent the Pauli equation 
for two-component wave functions and particles moving on a line or 
circle. Then n = k,s, with k = Xk/a labeling position (cf. (fTTTD . and 
s labeling the spin state. 

Stochastic trajectories for the discrete values of the system state n 
can be computed as follows. In a small time interval dt, state m will 
jump to n with probability P(m —>• n), where 

P(m n) = T nm dt. (14) 

The transition probability T nm is defined as 

Tnm — max(0, Jnm/Pm)i (15) 


where J is the probability current, 

2 

Jnm — —I ’n\{^ n Hnm'4 ) m) i 


(16) 


which is anti-symmetric in n and m, and P m is the usual probability 
for state m, 

Pm = CVW (17) 

The matrix H mn is the Hamiltonian in the n representation. Fur¬ 
ther details are in refs. [8] and [3][Ch.l9]. The above prescription for 
n-trajectories has the form of a guidance relation similar to ([3]) for 
(continuous) particle position: to compute the transition probabilities 
for time t, the wave function ijj n at time t must be computed from the 
wave equation (fT3l) . 

To obtain a wave function-free dynamics, the time-dependence of 
T nm should be computed without recourse to a guiding wave function. 
As in formulation F-III, the initial state should then be defined as an 
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ensemble of n- values, with prescribed probabilities P n . The stochastic 
trajectories along with the evolving probabilities P n (t) can then be 
computed as, 

P n {t) = P n (t ) ^2{T nm (t)P m (t) - T mn (t)P m (t )) (18) 


provided that the time-dependence of T nm can be computed from a 
relation like 

T nm (t) = F(T,P). (19) 

The function F may depend on the transition rates T and probabilities 
P, but must not refer to a pre-computed wave function as in ([15I16D . 

To discover the evolution equation (1191) . it is convenient to first 
find an evolution equation for the currents J nm , and then use (1T51) to 
convert it into an evolution equation for the transition rates T nm . The 
first step then is to rewrite the probability current (1 161) using a wave 
function in polar form (|2|), 


Jnm = -URnH n mRme-< Sn - S ^ /h ~ R m H rrm R n e i ^~ S ^l n ) ( 20 ) 

n 

and use this relation to solve for the wave function phase, or rather 
for 

0 nm = e -i(Sn-Sm)/h' ( 21 ) 

Since (1201) leads to a quadratic equation for 9 nm , there will be two 
solutions for a given J nm . However, consistency requires that 9 nm = 1 
when S n = S m , which implies that the solution takes the form, 


@nm — 


T ihJ n 

‘2Rm.Rn.Hn.rn. 


( 22 ) 


where 

[Hum] = + S ign (Re (H nm ))\H nm \ (23) 

and 

a nm = (4 RlRl - (24) 

It can be readily shown that the term in the square root is non¬ 
negative, hence a is real valued, with a nm = a mn > 0. The definition 
(1231) shows how the sign ambiguity is fixed; note that [H mn \ = [H nm ] 
and [H nm ] = H nm for a real valued Hamiltonian matrix. 

As stated above, a wave function-free dynamics, requires an expres¬ 
sion for J, which may depend on J, R (or P) and the Hamiltonian, but 
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not on the phases S. In such a formulation, the phases S (or equiv¬ 
alently, the 9) are derived properties, which can be back-computed 
from the evolving values of J and R. It should be noted that the ini¬ 
tial values for T, or equivalently J, cannot be chosen freely, since J nm 
must obey (1201) in which the R n and S n are the independent degrees 
of freedono 

The sign-choice implemented above assumes that in the initial state 
the real part of 9 nm is larger than zero. If, during subsequent time 
evolution 

ot-nm 0, (25) 

it may be required to flip the sign in (1231) such that 9 nm has a contin¬ 
uous time dependence. 

Using the result (1221) the time-dependence of J nm can be computed 
from 


Jn 


as 


= ^2 Re (X! R kHknHi 

k\n 


D A* A* 
nm^mV mn u nk 


k\m 


n Hum H m k Rk 9 nrn @mk) 


(26) 


J nm. — 


2 Rl 


U 


k\n 


nm / j Jnk' ~T2 [Hn 

n 


](%nm ^~' / [Hnk]c^nk') {*Tl Tl). (27) 

k\n 


The notation k\n in the summations indicates that only k values for 
which Hkn = is not equal to zero are included in the sum. 

To turn eq. (1271) into an evolution equation for the transition rates 
it is convenient to introduce an generalized transition rate T, defined 
as 


T 


nm 


Jnm/Rm: 


(28) 


where T nm is allowed to be negative, and 


Rnrn — Jnm/Pm (Jnm / Pm) Pm / Pn 


(29) 


Using m an d the continuity equation 


Pm — ^ ) Jmki 

k\m 


(30) 


2 It is assumed here that there are enough non-zero off-diagonal matrix elements in the 
Hamiltonian, such that for each independent S n (except one, because the overall phase is 
irrelevant) there is at least one non-zero element H nm , and hence J nm - 
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the final result is obtained as 


Trim — 2'^' nrn ('^'k\n'Tkn ^2k\mTkm ) 

+ 2 h 1 ^™. [Hnm] a nm(P n Sfc| n [Hkn]cXkn ~ P m J2k\m[Pkm\(Xkm) 

(31) 

With this result, a wave function-free formulation for the stochastic 
n-trajectories has been obtained, which is similar to formulation F-III. 
Specifically, the procedure to compute stochastic trajectories for the 
n-values is as follows. Define initial conditions for the trajectories of 
all n, by choosing probabilities P n and transition rates T nm . Since not 
all T nm can be chosen independently, it is actually easier to choose 

initial values for the phase S n and then use eqs. m - along with the 

1 /2 

chosen values for R n = P n - to compute the J mn and from these the 
T mn and T mn . In this way T mn will satisfy all required consistency 
conditions and subsequent evolution using (l32l) will preserve these 
conditions. Finally, a sufficiently small time step size dt must be 
chosen, such that also T mm dt = 1 — J2 n Pnmdt > 0 holds for all m. 
The time-dependence of the system is obtained by iteratively updating 
the probabilities P n and transition rates T nm , using 

P n (t + dt) = P n (t)+ J2 m (T nm (t)P m (t) 

Tmn C t)Pn(t )), 

T nm (t + dt) = T nm (t)+ T nm (t)dt, (32) 

T n m (t H - dt) — max(0, Tnmit H - Tt))^ 

where T is computed from eq. EH). The sign-choice implied in the 
[F/]-terms in (1311) should be the same as used for the previously time 
step, unless the cross-over condition (I25h is encountered, which might 
required to flip the sign in (1231) . 


4.1 Particle Moving on a Circle 

The update algorithm (|32l) is conceptually similar to formulation F-III 
for spin-zero particles in continuous (ID) space. To show this more 
explicitly, it is instructive to compute the update rules in the contin¬ 
uum limit for a particle moving on a circle, for which the Hamiltonian 
can be written as (cf. [8]), 

h 2 

)dm,n — M 2 (<5m+l,n + 5 mjn _|_i), 



Hmn — (Fn T 


Ma 2 
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with V n the potential, M the particle mass and a the grid spacing. 
In this case, the evolution equations of formulation F-III should be 
recovered. 

For a state in which the “particle” moves in the positive n direction, 
it is sufficient to compute the transitions rates from n to n + 1. This 
rate multiplied with the grid cell spacing a, can be interpreted as the 
average particle velocity, v n (cf. ([T2j0 . Using this change in notation, 
the generalized transition rates can be written as 


Tn+l,n — V n /a, Tn,n+l — V n P n /aP n + l- (34) 

Using this definition of velocity, the evolution equation (1311) turns into 

v n = 77 - v n (-v n -i + v n + v n - v n+ i) + (H-terms), (35) 

2a P n P n +i 

where “H-terms” indicate the second group of terms in (1311) , which will 
be considered shortly. Since the particle moves in one space dimen¬ 
sion, the summations in (1311) only have two terms, which are explicitly 
shown in (1351) . To investigate the continuum limit, a —>• 0, the veloci¬ 
ties are expanded to order a, v n ±\ = v n ± av' n and similarly for P n ±\. 
Substitution in (|35l) shows that 0(1) and 0(a) terms cancel, such that 
only the H-terms remain. 

Eq. (1331) shows that H n+ \ n = 0(l/a 2 ), from which it follows that 
a n +i,n = ‘ZRn+iRn + 0(a 2 ) and hence, 

Vn = ^ [Hn l 1,n] (Rn E ^M-^E^W) + 0(“)' 

n Un k\n+l k\n 

(36) 

Since the Hamiltonian matrix elements are real valued, the square 
brackets can be dropped and (1331) implies that 

[H n+ i, n ] = [H n , n+ 1 ] = - h 2 /2Ma 2 , [. H n>n \ =V n + K 2 /Ma 2 . (37) 


After substituting this in (f36l) and expanding the shifted terms R n +k 
to third order in a as R. n +k = Rn + akR' n + a 2 k 2 R” l /2 + a 3 k 3 R'" /6, 
one finds that 0(1), 0(a) and 0(a 2 ) terms cancel and the expected 
result emerges, 


Vn 




rfK 

2 MR n 


y + 0(0). 


(38) 
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4.2 Particle with Spin 1/2 

To keep this example simple, the particle is assumed to have a very 
large mass such that its motion can be ignored, and to only interact 
with a magnetic field B. Hence, the 2x2 Hamiltonian takes the simple 
form, 

H = ficrB = ^Ba x , a x = ^ J j , (39) 

where the x axis has been aligned with the magnetic field. Since 
n takes only two values, there is only one independent phase (the 
arbitrary global phase can be used to set one of the phases equal to 
0) and one independent probability (since P\ = 1 — P 2 ). Therefore 
there is also only one transition rate T i 2 = — T 21 P 1 /P 2 that has to be 
solved. Eq. (1311) therefore simplifies to 


Tv2 


5^< 1 + t ) + 


c^UHn? 1 
2 h 2 P 2 { Pi 



(40) 


Using [# 12 ] = fJ.B and 

ai 2 = (4 PM - ^TuP ?)'/ 2 


(41) 


it follows that 


T 12 = t\ 2 + 


2 n 2 B 2 


(2 



(42) 


Differentiating the continuity equation P 2 = —T 12 P 2 and eliminating 
T 12 using (fi 2 |) gives 


9.. 2 d2 

^ = -2-L(p 2 _i/ 2 ). 

This equation is easily solved as 


(43) 


P 2 (£) = (1 + cos( 27 f + 5 ))/2 = cos 2 ( 7 1 + <5), (44) 

with 7 = 2fj?B 2 /fi 2 and 5 a free integration constant. The generalized 
transition rate then follows as 


Ti 2 = 2 7 tan( 7 f + J). (45) 

It is straightforward to check that these results exactly agree with the 
results obtained from the wave function that solves the Schrodinger 
equation for this system. 
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5 Revisiting Entangled Spin States 

With the wave function free formulation as developed above, it is 
clear that also quantum non-locality and entangled spin states can 
be properly captured. This is the case, because there is a one-to- 
one relation between the dynamics of the probability held P n and 
(generalized) transition probabilities T nrn and the modulus R n and 
wave function phases S n - provided that the initial transition rates (or 
equivalently currents J nm are chosen in accordance with the constraint 
(1201) . Hence, the description of the measurement process, including 
non-local effects due to entangled quantum states that was presented 
in [1J [Ch.6,7] remains valid for the alternative formulation presented 
here. In formulation F-I the guiding wave function contains the extra 
spin state information, in the stochastic wave function-free dynamics, 
the state itself, which includes spin and position values as well as 
velocities and transition rates, contains all required information. Even 
though the evolution equations (|32l) may seem to be local, they of 
course describe the dynamics of a single point in the high-dimensional 
state space. An ensemble of such states can capture all required [3] 
correlations between the individual components of this state vector, 
also when these components refer to particle locations and spins at 
widely separated locations. 


6 Discussion 

The results of this paper extend previously developed formulations of 
quantum dynamics without wave functions to (non-relativistic) par¬ 
ticles with spin. Specifically the formulation developed in section [3] 
shows that the stochastic trajectories of the alternative formulation 
of quantum mechanics in terms of discrete “beables” [8], can be gen¬ 
erated without recourse to a guiding wave function. Exploring such 
an alternative might offer new insights. In particular it simplifies the 
“beable” world, or ontology, compared to a world in the Bohm-Bell 
formulation, F-I, which needs two different entities, particle locations 
with classical-like trajectories and a complex valued wave function 
that evolves according to the Schrodinger equation. Also, having an 
alternative set of evolution equations could lead to different, possi¬ 
bly more efficient numerical methods to solve and depict quantum 
dynamics. 

However, the stochastic dynamics developed above also has less 
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attractive features. Since the dynamics is computed in terms of the 
over-complete set of generalized transition rates, instead of the wave 
function phases, the initial conditions are constrained in a rather com¬ 
plex and unnatural fashion. Something similar is present in the wave 
function-free formulations F-II and F-III: here the particle trajectories 
have a constraint on the initial velocities (two or more dimensions). 
Even formulation F-I has a constraint on the initial probability density 
for particle positions, which must be chosen to be equal to the modulus 
squared of the wave function. However, here Bohm and others mm 
argue that these constrained initial conditions may not be required, 
since deviations will dynamically average out. In the current formula¬ 
tion, which lacks a guiding wave function, such a dynamical recovery 
of the interdependency of probability current and wave function phase 
seems very unlikely. 

Also the gain in simplicity compared to the Bohm-Bell formulation 
F-I is perhaps fairly meager: As in formulation F-III, one still needs to 
solve for the probabilities and transition rates in the full state space, 
i.e, an equally imposing task as solving the guiding wave function in 

F-I. 

Finally, there is the question if it is natural to adopt an integer spin 
value as an element of reality. In itself this is attractive, since it avoids 
having to argue that some quantum state properties are fundamen¬ 
tally different from others. I.e., one needs not argue that “location” 
is an intrinsic property, also in the microscopic world, but “spin” is 
not, and - unlike position - only emerges when macroscopic measure¬ 
ments are performed. Of course, as discussed in ref. [SJ, if spin would 
be granted a “beable” status, it is still not clear which spin: In the 
formulation developed above, eigenvalues of any spin operator can be 
“beables”. Choosing for example eigenvalues of cr z , the probabilities 
and transition rates are directly relevant for the evolution of the spin 
in this z-direction. Measurement result of a z can then be directly un¬ 
derstood, i.e., “read off” from the quantum mechanical state - as is the 
case for properties in a classical state. However, results of measure¬ 
ments of other (non-commuting) spins can no longer be understood 
directly and one has to invoke the same arguments as used in F-I 
([l][Ch.6,7] to understand how measurements of these spins will still 
agree with the results of normal quantum mechanics. 

It may be interesting to note at this point that one could resort 
to the view proposed in ref. [8j, that in fact all observables have a 
“beable” status. For spin this would imply that the particle has many 
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coexisting spin properties, for all possible spin operators (i.e., a dense, 
but still finite subset of all possible Hermitian and unitary 2x2 matri¬ 
ces). Quantum entanglement then plays out in a perhaps surprising, 
but very compelling way: In the typical set-up of two free moving par¬ 
ticles with entangled spins, the spin part of the Hamiltonian is zero, 
and hence the spin state remains unchanged while the particles move 
away from each other. Therefore it may seem that this system has the 
same (non)locality as a classical system in which some properties are 
correlated over long distance, simply because they started out that 
way. In order then to understand how the correct quantum results 
are obtained when performing spin measurements in arbitrary rela¬ 
tive directions - without evoking the decohering effect of a measuring 
apparatus - one must assume that the specific set-up of the measuring 
apparatus, specifically reads-off the property from the trajectory with 
a spin direction that aligns with the device. However, this cannot 
be done in a local fashion: the specific two particles spin realization 
is picked-up for which each of the two spin operators aligns with the 
corresponding measuring device. From a quantum trajectory point of 
view, the combined measuring apparatus reads-off a specific state from 
one trajectory, but this trajectory of course describes the correlated 
spin values for two widely separated particles. 

On the one hand, this coexisting states interpretation of quantum 
mechanics compellingly shows the non-local nature of “The Undivided 
Universe”; on the other hand - to quote Bell - “It seems an eccentric 
way to make a world”. 
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